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ABSTRACT 

We study the behavior of fracture in disordered systems close to the breakdown point. We 
simulate numerically both scalar (resistor network) and vectorial (spring network) models 
with threshold disorder, driven at constant current and stress rate respectively. We analyze 
the scaling of the susceptibility and the cluster size close to the breakdown. We observe 
avalanche behavior and clustering of the cracks. We find that the scaling exponents are 
consistent with those found close to a mean-field spinodal and present analogies between 
the coalescence of microfractures and the coalescence of droplets in a metastable magnetic 
system. Finally, we discuss different experimental conditions and some possible theoretical 
interpretations of the results. 

INTRODUCTION 

The breakdown of solids under external forces is a longstanding problem, that has practical 
and theoretical relevance. The way a material breaks, under the effect of an external electric 
field or under mechanical stress are closely related problems, due to the formal similarities in 
the underlying laws governing those phenomena. The first theoretical approach to fracture 
mechanics dates back to the twenties with the work of Griffith [l[ , who formulated a theory 
of crack formation, which is similar to the classical theory of nucleation in first-order phase 
transitions. Cracks grow or heal, depending on whether the external stress prevail over the 
resistance at surface of the crack. Similarly in bubble nucleation |2j , a critical droplet will 
form when the change in free energy due to the bulk exceeds that of the surface. Griffith 
theory assume the presence a single microcrack of a particular shape surrounded by an 
homogeneous medium, and therefore is not appropriate in disordered systems, where cracks 
can start from different positions and coalescence may take place. 

Spinodal nucleation ||, contrary to classical nucleation, is characterized by scaling prop- 
erties and fractal droplets. The spinodal point in fact has some characteristics of a critical 
point in second order phase transitions. The similarity between a solid driven to the thresh- 
old of mechanical instability and spinodal nucleation has been discussed in the past. Rundle 
and Klein Q], using a Landau-Ginzburg analysis of a single crack, showed that the crack 
growth obeyed scaling laws expected for spinodal nucleation. Selinger et al. || have shown 
by numerical simulations and mean-field theory of thermally activated fracture that the 
breakdown has the characteristics of a spinodal point. 

In this paper we concentrate on disordered media and we disregard the effect thermal 
fluctuations. The system is driven by an increasing external load to the point of global 
failure. It has been experimentally observed that the response, detected by acoustic emission 
(AE) measurements, to an increasing external stress takes place in bursts or avalanches 
distributed over a wide range of scale. Examples of this are found in foam glasses 0, fiber 
matrix composites H], concretes M, hydrogen precipitation || and volcanic rocks [IU|. We 
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observe a similar behavior for two dimensional discrete models. We show that the scaling 
behavior close to the breakdown is in quantitative agreement with the mean-field theory and 
it is suggestive of a first-order transition. The values of the scaling exponents are consistent 
with those found close to a spinodal point in thermally driven homogeneous systems. 

THE MODELS 



We study here two models, the random fuse model [if for electric breakdown and a spring 



network model Jl2| for fracture. In the random fuse model []TT[ each bond of a two dimen- 



sional lattice is occupied by a fuse of conductivity cr = 1, which burns when the current 
flowing in it exceeds a quenched random threshold. We consider a rotated square lattice 
with periodic boundary conditions in one direction. We impose a constant external current 
on the two other edges of the lattice. The currents in each bond are computed solving the 
kirchhoff equations. This step corresponds to the minimization of the total energy dissipated 
in the lattice 

£(M) = ~$>(A^, (1) 

i 

where (AV)j is the voltage drop in the bond i. We employ a multigrid relaxation algorithm 
with precision e = 1CU 10 . When all the currents are below the threshold we increase the 
current until the next bond reaches the threshold. The process is continued until a path 
of broken bonds spans the lattice and no current flows anymore. We chose a uniform 
distribution of thresholds, D 6 [0,2]. 



The second model is an elastic network [12 which has central and rotationally invariant 



bond-bending forces. The potential energy is 

£ = ^£(<NV i + ~ E(<%) 2 ^- (2) 

Z i Z <i,j> 

where 5r t is the change in the length of the bond i and 56ij is the change in the angle 
between neighboring bonds i and j. The constant <7j is equal to one if the bond is present 
and it is zero otherwise. A slowly increasing external stress is applied on all the edges and 
the lattice dynamics is obtained by numerically solving the equations of motion for each 
spring. Bonds break when stretched beyond a randomly chosen threshold. 

SIMULATION RESULTS 

The response of the model to the increase of the external force takes place in widely dis- 
tributed avalanches. The average size of the avalanches (i.e. the number of broken bonds) 



increases when the global failure is approached. We were able to show [13) using mean-field 



theory that the average avalanche size (m) diverges at the breakdown as 

(m) ~ (f c - /p 7 = 1/2- (3) 

where / is the stress or the current per unit length imposed on the lattice. We note that 
the same scaling law is expected close to a spinodal point, in the case of thermally driven 
first-order transitions. The macroscopic quantities of the system (i.e. elasticity) have a 
finite jump at the breakdown, indicative of a first-order transition. 
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Figure 1: a) The average avalanche size (m) is plotted as a function of / = I / L, the fit is 
done using the mean-field value 7 = 1/2. b) The "susceptibility" of the spring network with 
the mean field fit (7 = 1/2). The parameter <p is the fraction of bond that are not broken. 
The average avalanche size (m) is proportional to d<p/df. 

We confirm the validity of mean-field scaling by computer simulations of two dimensional 
models. For both models mean-field theory is obeyed remarkably well (see Fig. |I]a and 
Fig. |b). 

The reason for the observed mean-field behavior is probably due to the long-range nature 
of elastic interactions. The formation of cracks in those models takes place by the coalescence 
of several microcracks. This is confirmed by the behavior of the average crack size which 
does not diverge at the breakdown (see Fig 0). 

DISCUSSION AND CONCLUSIONS 

The breakdown of driven disordered media is described by scaling law which are reminiscent 
of those found close to a spinodal point. It appears that the behavior of a driven disordered 
system is similar to that of a thermally driven homogeneous system. This analogy is not 
too strict since the concept of metastability and spinodal are not well defined in the first 
case. 

Despite several experimental investigations of avalanche dynamics in fractures || [7]. [| 
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Figure 2: The average crack size for the fuse model as a function of the current for different 
systems sizes. The crack size does not seem to diverge at the breakdown. 

|], pT0[] , there is not a clear theoretical interpretation of the results. We believe that different 
experimental conditions can all give rise to similar scaling behavior, but the underlying 
physical mechanisms could be quite different. We can distinguish the following experimental 
setups: 

1. A solid driven by a constant stress rate can be described in the framework discussed 
in this paper. The system responds to the increase of the external load by AE bursts 
of increasing size 0, diverging at the point of global failure. It would be interesting 
to check if the scaling exponents agree with the mean-field theory. 

2. A solid subject to a constant load breaks because of thermal fluctuations. The AE 
is due to the formation of "droplets" and should be power law distributed close to 
the limit of stability (spinodal). Scaling exponents consistent with those of spinodal 
nucleation were observed in a recent experiment on cellular glass ||. To confirm this 
interpretation it would be necessary to study the scaling for different values of the 
applied load. 

3. A solid in a perfectly plastic state could respond to the increase of the external strain 
by a stationary AE signal. In this case one can interpret the results as a manifestation 
of self-organized criticality. Such a behavior was shown in numerical models [0], |T5|j , 
but to our knowledge it has not yet been observed in experiments. 

We believe that extensive and systematic experiments along these lines can resolve these 
longstanding problems. 
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